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We propose a technique for performing quantum state tomography of photonic polarization- 
encoded multi-qubit states. Our method uses a single rotating wave plate, a polarizing beam splitter 
and two photon-counting detectors per photon mode. As the wave plate rotates, the photon counters 
measure a pseudo-continuous signal which is then Fourier transformed. The density matrix of the 
state is reconstructed using the relationship between the Fourier coefficients of the signal and the 
Stokes' parameters that represent the state. The experimental complexity, i.e. different wave plate 
rotation frequencies, scales linearly with the number of qubits. 
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Quantum state preparation is an essential ingredient 
in the realization of quantum technologies such as quan- 
tum computing p], quantum cryptography p] and other 
quantum information protocols [3]. A crucial aspect of 
reliable state preparation is the ability to accurately char- 
acterize the state of a quantum system. To this end, 
quantum state tomography (QST) allows the reconstruc- 
tion of a state's density matrix from measurement statis- 
tics accumulated through repeated independent measure- 
ments of multiple identically-prepared systems |3H6] . 

In linear-optics, where quantum information is en- 
coded in the polarization of a single photon, different 
measurement settings are realized with a combination of 
linear optical elements such as wave plates, beam split- 
ters and polarizing beam splitters, followed by photon 
counting. QST was first accomplished in such systems 
by White et al. [7J, where the measurement settings 
corresponded directly to the Stokes' parameters used to 
characterize the polarization state of the classical elec- 
tromagnetic field |S]. Later it was suggested that an 
over-complete symmetric six-measurement set [S] or an 
informationally-complete symmetric four-measurement 
set [9lfT3] be used for improved performance. Other ex- 
tensions, such as those considering optimal experimental 
design under realistic technical constraints [2J [15], or 
modifications due to inaccessible information [T&H22] or 
preferable measurements choices [SJ [TTJ [T^l have 
also been considered. 

To date, all implementations of QST of photonic 
polarization-encoded qubits have utilized either multiple 
wave-plates and/or multiple beam splitters per qubit. We 
propose a technique that uses only one wave plate and 
one polarizing beam splitter (PBS) per qubit mode. Each 
mode m is incident on a single wave plate rotating at fre- 
quency f2 m followed by a polarizing beam splitter (PBS). 
Photon counters at the output ports of the PBS measure 
a pseudo- continuous signal and the state is reconstructed 
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from the Fourier coefficients of this signal. The experi- 
mental complexity of this method scales linearly with the 
number of qubits in terms of the number of settings re- 
quired (i.e. wave plate rotation frequencies) rather than 
exponentially, as is the case with QST that uses discrete 
measurement settings. Similar techniques that rely on 
rotating wave-plates are used in classical optics to deter- 
mine the polarization state of the electromagnetic field 
[31] , In the context of non-classical light, Fourier spec- 
troscopy has been used to characterize the joint spectrum 
of photons [32] . 

The remainder of this paper is organised as follows. 
In Section [ij we give a brief review of QST of multi- 
qubit states. In Section[nJ we introducing our scheme for 
Fourier transform tomography (FTT) and provide exam- 
ples for one and two qubits. In Section |III| we provide 
concluding remarks 

I. QUANTUM STATE TOMOGRAPHY 

Tomography is the process of constructing a represen- 
tation of an object by imaging it in different sections. In 
quantum state tomography, we aim to construct a rep- 
resentation of a quantum state p from different measure- 
ment outcomes. An n-qubit system is specified by 4™ — 1 
real parameters. We therefore require at least this many 
outcomes of linearly independent measurements to spec- 
ify P- 

The probability of obtaining measurement outcome j, 
given a measurement operator Mj, is given by 

Pj = (Mj) = Tr[pMj] = J-, (1) 

where rij is the number of counts and Afj is a constant 
dependent on the detector efficiency and duration of data 
collection. In a polarization-encoded linear optical sys- 
tem, any projective measurement can be realized with 
a quarter-wave plate, a half-wave plate and a polarizing 
beam splitter, as shown in FIG. [l]a). A popular choice 
corresponds to the three Pauli operators. 
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FIG. 1: (Color online) Schematic diagrams of: a) typical 
QST set-up which uses a combination of quarter- and half- 
wave plates to perform arbitrary-basis measurements; b) FTT 
set-up which uses one rotating wave plate; c) multi-qubit FTT 
which uses one wave-plate per qubit mode m, rotating at 
frequency f2 m . 



We can always write the density matrix of an n-qubit 
system in terms of Hermitian operators &i 



P = 



1 3 
2™ ^ 

i 1: ...,i n —0 



(2) 



where &o = (H \ + \V) (V\ is the identity operator and 
<xi_3 are the Pauli operators: <j\ = \H)(V\ + \V)(H\, 
& 2 = i{\V){H\ - \H){V\) and a 3 = \H) (H\ - \V)(V\. 
The coefficients Si lt ,,. t i n = Tr{p(<Ji 1 ® ••• ® <7,- n )] com- 
pletely characterize the state. j n are normalized 
generalizations of the classical parameters introduced by 
Stokes in 1852 [3J, and will hereafter be simply referred 
to as Stokes ' parameters. 

Combining Equations ([T]) and ^ , we find a linear rela- 
tionship between Stokes' parameters and the probability 



1 

2" 



E a*. 

il,...,i n = 



®& in Mj\, (3) 



where Mj acts on the entire multi-qubit system. By mak- 
ing 4™ — 1 linearly independent measurements, it is pos- 
sible to solve for Stokes' parameters and reconstruct the 
density matrix according to Equation This can be 
achieved through a variety of methods, including simple 
linear inversion, least-squares estimation or the popular 
maximum likelihood estimation method 33J . Alterna- 
tively, one can look to a growing number of exciting new 
techniques such as the forced purity routine [53], Bay- 
sean mean estimation [35], compressed sensing [3Jj], von 
Neumann entropy maximization [37], hedged maximum 
likelihood estimation [38], minimax estimation [39], and 
techniques that focus on reconstructing the state with 
reliable error bars [3H| and confidence regions [4"T] . 



II. FOURIER TRANSFORM TOMOGRAPHY 

In this section, we show how the quantum state of a 
multi-qubit system can be represented by a single joint- 
probability signal and how the measurement of this signal 
enables the reconstruction of the quantum state. 

In our proposal, identical copies of the state are pre- 
pared and subsequently pass through a series of optical 
elements. For a multi-photon state, each photon mode m 
is incident on a single wave plate rotating at frequency 
O m followed by a polarizing beam splitter (PBS). Pho- 
ton counters at the output ports of the PBS continuously 
measure the intensity, which can be processed to recover 
Stokes' parameters. A schematic of this setup is shown 
in in FIG [I] b) for a single qubit and [I] c) for multiple 
qubits. For multiple qubits, the signal measured is a 
"coincidence intensity" corresponding to the joint prob- 
ability of detecting photons at each PBS. 

The time-dependent single-qubit projection- valued 
measure (PVM) associated with the probability of detect- 
ing a photon in the horizontal or vertical output modes 
of each PBS is given by {M%(t),M%{t)} where 

M^(t) = Ul\a)(a\U m , (4) 
for a = H,V , where m labels the qubit mode and 

U m (t) = cos f^j & - ism P^J v m (t) ■ a (5) 

is the unitary operator associated with a wave plate in 
mode m. U m (t) rotates the operators \a){a\ on the Bloch 
sphere by an angle /?, about the vector 



v m {t) = cos(ui m t)k + sin(w m t)j , 



(6) 



where k and i are unit vectors in Euclidian space (defined 
by the axes in FIG[2| and v ■ a = v\<j\ + V2&2 + As 
the wave plate rotates about the beam-axis at frequency 
Q m in real space, v m (t) rotates about the y-axis in Eu- 
clidian space at frequency uj m = 2Q m . We assume that 
the fast axis of the wave plate is aligned at degrees to 
the horizontal as defined by the polarization of the pho- 
tons. A phase factor can be included in v m (t) to account 
for different initial alignment of the wave plate. The re- 
sulting projector M^f (t) traces out a figure-8 path on the 
Bloch sphere, as shown in FIG [2] The retardance of the 
wave plate determines the size of the figure-8. 

To characterize an n-qubit state, one measures a joint 
probability of detecting a photon in the H mode of each 
PBS. This is given by 



>(*) = 4 E Sn. 



.Xl.ii • ■ • Xr, 



ii,...,i n =0 



where 



Tr[&iM%(t)] 



(7) 



(8) 
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FIG. 2: (Color online) Path traced out by M^\t), defined 
in Equation Q, for: /3 = tt/4 (green, dotted); /3 = ty/2 (blue, 
dashed); and /3 = ll7r/15 (orange, solid). 



and therefore 



Xm,0 — 1 

Xm.i = s 2 sin (2oj m t) 

Xm,2 = 2cSSU\(uj m t) 

Xm,3 = c 2 + s 2 cos (2w m t) 



(9a) 
(9b) 
(9c) 
(9d) 



where c = cos(/3/2) and s = sin (/?/2). 

Note that the choice of analyzing the signal from mode 
H rather than mode V is arbitrary and typically both 
modes will need to be measured to ensure normalised 




FIG. 3: (Color online) The equally weighted variance (EWV) 
|42| assesses the noise immunity of the wave-plate. Here we 
plot EWV x N for one period as a function of /3 (orange, solid) , 
where iV is the total number of bins per period. A smaller 
EWV is associated with better immunity to noise. The best 
noise immunity occurs at j3 ~ 2.27 ~ ll7r/15 for EWV « 
41.7/N. The black, dotted line shows EWV opt = 40/N, at- 
tainable by optimal tomographic schemes such as those that 
measure the Pauli matrices or the SIC-POVM [TTJ1 fTTl ITS]. 
The blue, dashed line shows the EWV for a QWP Q3 = tt/2), 
EWV QWP = 84/iV. Note that the EWV is dimensionless. 



probabilities. 

Without loss of generality, we restrict < cji < • • • < 
Lu n . For two qubits, u; 2 = rwi where r > 1. If 7' is 
an irrational number, the signal does not have a finite 
period. If r is a rational number, we can write r = p/q, 
where p and q are integers. In this case, the period of the 
two-qubit signal is given by 



T(r) = 



2irq 
wigcd(p, q) 



(10) 



where gcd(p, q) is the greatest common denominator of 
p and q. For n > 2, the period of the signal can be 
determined via recursion. A shorter period is favourable 
from an experimental perspective which, for a constant 
loi, occurs when r is an integer. The lowest integer that 
ensures sufficient Fourier coefficients to solve for Stokes' 
parameters is r = 5. 

In practice p n will not be a continuous function of time 
but rather a discretized approximation. The discretized 
signal will be divided into time bins, with M coincidence 
counts in each bin. The number of time bins per pe- 
riod, N, must be at least the Nyquist rate, i.e. twice the 
highest frequency contained within the signal, to avoid 
aliasing. 

The discrete probability in bin Tj will be given by 



nH...H 

N 



(11) 



where N — Ylki k 2 n=H v n ki....,k 2 ™ i n ^ s number of 
qubits and ...,fc 2 „ is the number of coincidence counts 

for a given projector M^ 1 (£) ® • • • ® M^ ! " (t) . 

In principle, j3 can take on any value other than an 
integer multiple of n. However in practice, some values 
will be more susceptible to noise than others. We use 
the equally weighted variance (EWV) [42] , which assesses 
the noise immunity of the wave-plate, to show that (3 
ll7r/15 is most immune to noise, as defined in [32] • A 
plot of the EWV is shown in FIG [3] Such a retardance 
can be achieved with an off-the-shelf wave plate designed 
for a wave length different to that of the experiment. 

In the remainder of this section, we provide specific 
examples for one- and two-qubit states. 



A. Example: one qubit 

For a single qubit, the signal is given by 
1 3 



i=0 

S Si 

— + —s 2 sin (2uit) + S2CS sin (oj\t) 
g 

+ y (c 2 + s 2 cos (2wii)) . 



(12) 
(13) 
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FIG. 4: The probability signal p n , as a function of revolutions of the first wave plate, and Fourier coefficients af and 6/ for the 
state: a) pi given in Equation (181; and b) \tp)2 given in Equation (231. The label / denotes the term in the Fourier series and 



corresponds to the subscript in Equations (1141) and (122 



This can be written as 



parameter d, such that 



Pi{t) = y + b i sin(wit) 

+ a 2 cos (2ujit) + b 2 sin (2u>it) 

where the Fourier coefficients are given by 

a = So + S 3 c 2 ; 
h = S 2 cs ; 



«2 



V: 



. Si 2 
b 2 = -js . 



(14) 

(15a) 
(15b) 

(15c) 
(15d) 



where c = cos(/3/2) and s = sin(/3/2). Linear inversion 
of Equations (15) gives the Stokes' parameters in terms 



of the Fourier coefficients: 



So — ao — 



Si 

s 2 
S3 



262 

cs 
2a 2 

c 2 



(16a) 
(16b) 
(16c) 
(16d) 



Substitution into Equation ([2]), gives the density matrix 
in terms of the Fourier coefficients: 



\ (a + 2a 2 ) 



Pi — \ 63 1 „• fei 1 /„ 1 o„ \ 2a 2 



+ *2^ \{a Q + 2a 2 )- 



(17) 



As an example, consider a single-qubit state that 
has experienced depolarizing noise, characterized by the 



pi = d|^)i<^|i + (1 - djd-.l^i^lid-* 



Specifically, let's consider 
1 



\H) 



-J7T/4 



1*0) 



(18) 



(19) 



and d = 0.1. A rctardance of /? = ll7r/15 produces the 
signal shown in FIG[4]a). Performing a fast Fourier trans- 
form (FFT) of the discretized signal yields the Fourier 
coefficients in FIG [1] a) . The coefficients a f and b f corre- 
spond to the real and imaginary parts of the list gener- 
ated by the FFT respectively. Inserting the coefficients, 
a Q = 1, bi = 0.210, a 2 = and b 2 = -0.236 into the 



density matrix in Equation (17) gives 



Pi = 



0.5 

-0.283 + 0.283i 



-0.283 - 0.283i . } 
0.5 I ' [ Uj 



which corresponds to the density operator in Equation 
p} for d = 0.1. 



B. Example: two qubits 

For two qubits, the joint probability of detecting a pho- 
ton in the horizontal output ports of each PBS is given 
by 



P2{t) = 2^ — 7-Xl,iiX2,i a 



(21) 

ii ,i2=0 

y + ( a f cos ( w /*) + b f sin(wy*)) , (22) 
f=i 
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where in the second line, we have written the signal in 
terms of its Fourier coefficients. The extent of the sum- 
mation depends on the specific choice of relative frequen- 
cies, and uj'j are functions of uj\ and cj 2 from the set 
{wi, u>2, 2cji , 2w 2 , wi ±u>2, 2w 2 , 2wi ±uj 2 , 2lox ± 2oj 2 }- 
The elements of this set are not necessarily in order of 
size, and to relate them to uj'j one needs to consider ex- 
plicit values for lo\ and ui2- 

As an example, consider the two-qubit state 

\^) 2 = ^(\H)\V) + \R)\L)). (23) 



Wave plate retardances of (3 = ll7r/15 and a frequency 
ratio r = W2/W1 = 5 produces the signal shown in FIG 
[1] The explicit expression for P2(t) for this choice of 
measurement settings can be found in Appendix [X] Per- 
forming a fast Fourier transform (FFT) of the discretized 
signal yields the Fourier coefficients in FIG [4] b) . 

Linear inversion of the expressions for the Fourier coef- 
ficients in terms of the Stokes' parameters, given in Equa- 
tions ( |A3[ ), followed by substitution into Equation ([2|, 
along with the Fourier coefficients determined from the 

nalin FiG[|b), gives the reconstructed density matrix 



/ 0.125 0.25 + 0.125i — 0.125z 0.125 

0.25 - 0.125i 0.625 -0.125 - 0.25i 0.25 - 0.125i 

0.125i -0.125 + 0.25z 0.125 0.125z 

V 0.125 0.25 + 0.125z —0.125? 0.125 




which corresponds to the density operator p 2 — \ip)2 (tp \ 2 
where \ip)2 is defined in Equation (23). In general, given 
separable qubits, P2 factorizes into a product of p\ for 
each qubit. 



(USRA) and the DARPA (QuBE) program. 



III. SUMMARY & CONCLUDING REMARKS 

We presented a scheme for performing quantum state 
tomography of photonic polarization-encoded multi- 
qubit states. The scheme is simpler than standard to- 
mographic protocols in that only one wave plate and one 
polarizing beam splitter is required per photon mode. 

In this scheme, photon-counting detectors measure a 
pseudo-continuous time-dependent joint probability as 
the wave plates rotate at frequency fl m . The Fourier co- 
efficients of the signal give the Stokes' parameters which 
describe the state. For a single qubit, the optimal wave 
plate retardance is /? ~ ll7r/15. 

This technique reduces the number of required optical 
elements and the experimental complexity scales linearly 
with the number of qubits, in terms of the number of 
settings required (wave plate rotation frequencies) rather 
than exponentially, as is the case with QST that uses 
discrete measurement settings. 

An open question is whether the representation of a 
quantum state as a continuous signal will provide intu- 
itive means for establishing certain properties of the state 
such as its entanglement. 
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Appendix A: Probability signal and Fourier 
coefficients for two-qubit state 

In this appendix, we give the probability signal for 
the specific two-qubit example described in Section [IIB| 
We also provide expressions for the Fourier coefficients in 
terms of the Stokes' parameters, as well as the inverted 
expressions for the Stokes' parameters in terms of the 
Fourier coefficients. 

The signal probability for the state 

\1>) 2 = ±(\H)\V)+\R)\L)), (Al) 

with a frequency ratio r = ^2/^1 = 5, is 

P2(t) = + bi sin (uJit) + b 2 sin (2uit) + 63 sin (3wii) 

+ 65 sin (5uiit) + bj sin (Juiit) + bg sin (8u>it) 
+ bg sin (9uJit) + b\o sin (lOcjii) + bu sin (llwit) 
+ 612 sin (12wit) + a 2 cos (2ujit) + 03 cos (3wii) 
+ a 8 cos (8uJit) + a 4 cos (4^!^) + a 6 cos (6uit) 
+ 0,7 cos (7wii) + ag cos (9u>it) + aio cos (lOwit) 
+ an cos (llojit) + ai2 cos (12wi£) , 

(A2) 



G 



where the Fourier coefficients are 



(c 2 (c 2 5 3 , 3 + S0.3 + Ss,o) + 5*o,o) /2 
s 2 (c 2 5 3 , 3 + S 3 ,o) /4 
-a 7 = cs 3 5i, 2 /4 
cV5 2 , 2 /2 



no 



0*3 
«4 

aio 

012 = s 4 (5*3,3 - 5 M ) /8, 



= -an = cs 3 5 2 , i/4 
s 2 (c 2 5 3 , 3 + 5 0>3 ) /4 



(A3a) 
(A3b) 
(A3c) 
(A3d) 
(A3e) 
(A3f) 
(A3g) 
(A3h) 



and 





= cs 


(c 2 5*2,3 + 5 2 , ) /2 


(A4a) 


b 2 


= s 2 


(c 2 5i, 3 + 5 1)0 ) /2 


(A4b) 


h 


= b 7 


= cs 3 5 3 , 2 /4 


(A4c) 


h 


= cs 


(c 2 5 3 , 2 + S Q , 2 ) /2 


(A4d) 


b 8 


= s 4 


(53,i — 5i, 3 ) /8 


(A4e) 


b 9 




& u = -cs 3 5 2 , 3 /4 


(A4f) 


ho 


= s 2 


(c 2 5 3 ,i + 5 ,i) /4 


(A4g) 


&12 


= S 4 


(5 1 , 3 + 5 3 ,i)/8, 


(A4h) 



where c = cos(/3/2) and s = sin(/?/2). Inverting Equa- 



tions (A3) and ( |A4[ ) gives 

5 ,o = 4c 2 ((o 8 + aia) c 2 - (a 2 + 010) s 2 )/s 4 + 



2a n 







(A5a) 


5o,i 


= -4((& 8 + &i2)c 2 -6ios 2 )/ S 4 


(A5b) 


5o,2 


= 265/cs - 46 3 c/s 3 


(A5c) 


5o,3 


= 4ai /s 2 - 4 (a 8 + a 12 ) c 2 /s 4 


(A5d) 


5i,o 


-4((6 8 -6i 2 )c 2 + b 2 s 2 )/s 4 


(A5e) 


5i,i 


= 4(a 8 - ai 2 )/s 4 


(A5f) 


5l,2 


= 4a 3 /cs 3 


(A5g) 


5l,3 


= -4(& 8 -&i 2 )/s 4 


(A5h) 


5 2 ,o 


= 4& 9 c/s 3 + 26i/cs 


(A5i) 


52,1 


= 4a 9 /cs 3 


(A5j) 


52,2 


= 2a 4 /c 2 s 2 


(A5k) 


52,3 


= -4& 9 /cs 3 


(A51) 


5 3 ,o 


= 4a 2 /s 2 - 4 (a 8 + ai 2 ) c 2 /s 4 


(A5m) 


5*3,1 


= 4(& 8 + 6i 2 )/.s 4 


(A5n) 


53,2 


= 46 3 /cs 3 


(A5o) 


53,3 


= 4(a 8 + ai 2 )/s 4 . 


(A5 P ) 



Substituting the Stokes' parameters, along with the 
Fourier coefficients in FIG [4] b) , into Equatio n (|2| , gives 
the reconstructed density matrix in Equation ( |24[ ) which 
corresponds to the density operator p 2 = \ip) 2 (ip\ 2 where 
\tp} 2 is defined in Equation (Al). 
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